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ABSTRACT 



Context. Since most high- and intermediate-mass protostars are at great distance and form in clusters, high linear resolution observa- 
tions are needed to investigate their physical properties. 

Aims. To study the gas in the innermost region around the protostars in the proto-cluster IRAS 05358-1-3543, we observed the source 
in several transitions of methanol and other molecular species with the Plateau de Bure Interferometer and the Submillimeter Array, 
reaching a linear resolution of 1 100 AU. 

Methods. We determine the kinetic temperature of the gas around the protostars through an LVG and LTE analysis of their molecular 
emission; the column densities of CH3OH, CH3CN and SO2 are also derived. Constrains on the density of the gas are estimated for 
two of the protostellar cores. 

Results. We find that the dust condensations are in various evolutionary stages. The powerhouse of the cluster, mmla, harbours a hot 
core with T ~ 220 (75 < T < 330) K. A double-peaked profile is detected in several transitions toward mmla, and we found a velocity 
gradient along a linear structure which could be perpendicular to one of the outflows from the vicinity of mmla. Since the size of the 
double-peaked emission is less than 1 100 AU, we suggest that mmla might host a massive circumstellar disk. The other sources are 
in earlier stages of star formation. The least active source, mm3, could be a starless massive core, since it is cold (T < 20 K), with a 
large reservoir of accreting material (M ~ 19 Mq), but no molecular emission peaks on it. 

Key words. Stars: formation; stars: early type; stars: individual (IRAS 05358-1-3543); ISM: lines and bands; ISM: molecules 



1. Introduction 

The last decade has seen significant progress in the under- 
standing of how high-mass stars form. Large samples of mas- 
sive young stellar objects (YSOs) were studied with single- 
dish telescopes, to investigate their physical properties through 
the analys is of their (sub'Imm continuum and molecular emis- 
sion (e .g.. iMohnari et al.lll996l l 19981 120001; IWalsh et al.lll997 , 



19981 Il999t iHatchell et al.1 I2OOOI; IZhang et al.l 12001 



2005 



Sridharan et al.' '2002'; Beuther et all l2002bllcl; iFaundez et all 
2004; Wi lliams et al. 2004. 20051). However, an intrinsic feature 
of high-mass stars is that they form in clusters, and that most of 
them are at large (several kpc) distances. Therefore, single-dish 
studies, as valuable as they are, lack the necessary spatial res- 
olution to resolve single protostars and study the inner regions 
where high-mass star formation takes place. Interferometric ob- 
servations started shedding light into the complex nature of high- 
mass star forming regions with the adequate spatial resolution 
(e.g., Cesaroni et al. 1997, 1999, 2005; Wvrowski et al. 199^ 
Beuther et al.ll2002al 120051: IHunter et al.ll2006t IShepherd et all 
20001) . However, the number of massive YSOs studied at high 
resolution is still too small to establish the general properties of 
the dense cores where massive stars form on a statistical base. 

In this paper, we present an interferometric analysis of the 
high-mass star forming region IRAS 05358+3543 at (sub)mm 
wavelengths in several molecular transitions. IRAS 05358-1-3543 



(also known in literature as S2331R) is part of a sample of 
69 high-mass protostellar objects s tudied in great detail in 
recent years (Sridharan et al. 2002; Beuther etal. 2002b' cS 
Williams et al. 200 4, .200 5.; Fuller et al. 2005). At a distance of 
1.8 kpc iSnell et al.l ll990fc . IRAS 05358H-3543 has a bolometric 
luminosity of 6300 Lq; strong high- mass star formation activity 
is evidenced by maser emission (see lMentenl [T99 1 ; Tofani et aD 
'1995; Minier et al. 2000) and outflow activity (Snell et al.ll 19901 
[Beuther et al. 2002a,) . Previous interferometric observations by 
iBeuther et alJ (l2002ah resolved three dust condensations (mml, 
mm2 and mm3) within an area of 9" x4" (~ 17 100 x 7 200 AU), 
and revealed at least three outflows in CO and SiO, the most 
prominent of which is more than a parsec in length, and massive 
(M > 10 Mq). Two of the three identified outflows originate 
from the vicinity of mml, which is probably the main power- 
house in the region. 

To zoom in on the innermost region around the protostars, 
and study the physical properties of the individual potentially 
star-forming cores, we caiTied out a comprehensive program to 
observe the region at high spatial resolution with the Plateau 
de Bure Interferometer at 97 GHz and 241 GHz, and the 
Submillimeter ArrajO at 338 GHz. The new observations reach a 
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Table 1. Positions of the four dust condensations in 



Source 


R.A. [J2000] 


Dec. [J2000] 


mmla 


05:39:13.08 


35:45:51.3 


mmlb 


05:39:13.13 


35:45:50.8 


mm2a 


05:39:12.76 


35:45:51.3 


mm3 


05:39:12.50 


35:45:54.9 



resolution down to 0.6", corresponding t o ~ 1100 AU at the dis- 
tance of the source. iBeuther et al.l ( l2007ah studied the continuum 
emission of this dataset. They identified four compact protostel- 
lar sources in the region; mml is resolved into two continuum 
peaks, mmla and mmlb, with a projected linear separation of 
~ 1700 AU. A mid-infrared source (jLongmore et al. 2006), and 
a compact 3.6 cm continuum source (Be uther et alji2007al) co- 
incide with mmla, which is also associate d with the class II 
methanol masers detected by IMinier et aTl (l2000l) . The previ- 
ously identified source mm2 resolves into several sub-sources; 
however, on ly one of them (mm2 a, according to the nomencla- 
ture used bv 'Beuther e t al.ll2007al) is a protostellar source, while 
the others are probably caused by the outflows in the region. 
The third source mm3 remains a single compact core even at 
the highest spatial resolution. In Table [T] we report the posi- 
tio ns of the four sources identified in the continuum emission 
by iBeuflieretai] (l2007 a) . 

In this paper, we discuss the spectral line observations 
comple menting the continuum data discussed by iBeuther et al.l 
(l2007ah . In section ^ the different observations are presented. 
In section ^ we discuss our results, and analyse the ex- 
tended emission of low excitation molecular transitions ( ^3. It . 
as well as the molecular spectra at the positions of the dust 
condensations 03.2b . Finally, in section ^we derive the physi- 
cal parameters of the gas around the protostars from the analysis 
of their spectra. In the following sections, we use the term pro- 
tostar for young massive stellar objects which are still accreting 
material from the surroundings, independently whether they al- 
ready started burning hydrogen or not. 



2. Observations 

2.1. Plateau de Bure Interferometer (PdBI) 

IRAS 05358+3543 was observed with the IRAM Plateau de 
Bure Interferometer in two diff'erent frequency setups, in 2003 
and 2005. A first frequency setup was performed in four tracks 
between January and October 2003 in the BC and D configu- 
rations of the array. Two observations, in January and October, 
were performed with only 5 antennas; on January 12th, 2003, 
the array was configured in a special BC combination. The 
3 mm receivers were used in single-side band mode and tuned to 
96.6 GHz; the 1 mm receivers, in double-side band mode, were 
toned to 241.85 GHz (USB). At 3 mm, the C^'^S fine and the 
torsionally excited li^ —> CH3OH quartet were covered us- 
ing two correlator units of 80 MHz bandwidth. One 320 MHz 
unit was placed to obtain a continuum measurement at 3 mm. 
The 5i ^ 4i Vt = CH3OH band, the 2i,i ^ 2^2 HDO and 
the 5 — » 4 SO2 fines were observed with four units of 160 MHz 
bandwidth, which were also used to obtain a continuum mea- 
surement at 1.3mm. The same configuration of the correlator 
units allowed the observation of 13a ^"^k CH3CN band in 



the LSB at 239 GHz. The observations were performed in the 
MOSAIC mode, with seven fields covering the whole source in 
a hexagonal pattern with a centre field (see Fig.[T]). 

The second frequency setup was observed on February 4th 
and February 7th, 2005, in the A and B configurations. The 3 mm 
receivers were used in single-side band mode and tuned again 
to 96.6 GHz, with the same configuration of the correlator de- 
scribed before. The 1 mm receivers, in double side-band mode, 
were tuned to 241.2 GHz (USB) with three 160 MHz units cov- 
ering the 5i -> 4a Vt = 1 CH3OH band. 

Bandpass calibration was done with 0420-014, 3C454.2 and 
NRAO150 for the BCD configuration, with 3C84 for the AB 
data. NRAO150, 0528+134 and 3C273 were used as flux cal- 
ibrators of the BCD and AB data, respectively. Phase and am- 
plitude calibration was done via observations of 0528+134, 
0552+398, 0529+483 and J0418+380. Measured system tem- 
peratures were between 100 and 260 K in the 3.1 mm band for 
both setups. On February 4th, the system temperatures in the 
1.2 mm receivers ranged between 240 and 420 K, with one re- 
ceiver measuring system temperatures of 900 K. On February 
7th, and during the observations of the first frequency setup, the 
system temperatures in the 1.2 mm bands were higher, between 
400 and 1000 K, due to less favourable weather conditions. 

For both frequency setups, the phase centre was 
a20oo=05''39"'13 '.07, 520()o=+35°45'50".5, with vlsr=-17.6 km 
s '. Details on the spectral resolutions, the synthesised beam 
sizes and the main spectral lines per frequency units are given in 
Table |2] For the AB configuration the baselines range between 
30 and 400 m, for the BCD configuration between 20 and 
320 m. Therefore, at 1.2 mm any source structure larger than 6" 
is filtered out in the AB observations, and larger than 13" in the 
BCD data. 

The data calibration and the imaging were performed with 
the CLIC and MAPPING softwarefl For the PdBI and the SMA 
data, line free channels were averaged to produce continuum im- 
ages, which were then subtracted from the line data in the visi- 
bility planeQ 

2.2. Submillimeter Array 

We observed IRAS 05358+3543 with the SMA on November 
11th, 2004, at 348 GHz (865 jum) in the compact configuration 
with seven antennas, and on January 15th, 2005, in the extended 
configuration again with seven antennas in the array. However, 
the data from one antenna in the extended configuration were 
unusable reducing it to six effective antennas for that configura- 
tion. Due to problems at the correlator during the observations in 
the extended configuration, the frequency of the 7^6 C^'^S line 
was not covered. However, the 7 — > 6 C^'^S map was produced 
from the compact configuration data. The projected baselines 
ranged between 13 and 223 m. The short baseline cutoff implies 
that source structures > 16" are filtered out by the observations. 
The phase centre of the observations was a'2()oo=05''39'"13'.07 
and 52()00=+35°45'51".2 with a Vki- - -17.6kms"^ Bandpass 
calibration was done with Jupiter, Uranus, Callisto, and 3C279. 
We used Callisto and 3C279 for the flux calibration which is 
estimated to be accurate within 20%. Phase and amplitude cali- 
bration was done via frequent observations of the quasar 3C1 1 1 



Astronomy and Astrophysics and is funded by the Smithsonian 
Institution and the Academia Sinica. 



^ http://www.iram.fr/IRAMFR/GILDAS 

* The calibrated UV tables related to the PdBI and SMA observations 
are available in electronic form at the CDS. The integrated intensity 
map of the 5o ^ 4o v, = CH3OH-A line with the 30 m IRAM is also 
available in electronic form at the CDS. 
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Table 2. Observational parameters. 



Main line 




Centre freqtiency 


Configtiration 


rlrrJW 


r.A. 




r.m.s. 






(GHz) 




(") 


(°) 


(km s ) 


(Jy/beam) 


C'^'S 2 ^ 1 




96.44 


BCD 


4.21 X 3.07 


66 


0.5 


0.01 


CH3OH 2t It V, 


= 1 


96.51 


BCD 


4.21 X 3.06 


66 


0.5 


0.01 


CH3OH 2t ^ U V, 


= 1 


96.51 


AB 


1.85 X 1.36 


26 


1.0 


0.004 


CH3OH 2j. -» Ij. V, 


= 


96.72 


BCD 


4.19 X 3.04 


65 


8 


0.01 


CH3CN 13* ^ ilk 




238.90 


BCD 


2.12 X 1.26 


77 


0.8 


0.03 


C3*S 5^4 




241.08 


AB 


0.77 X 0.55 


13 


1.5 


0.01 


CHiOH 5 k 4k V, 


= 1 


241.30 


AB 


0.77 X 0.55 


13 


1.5 


0.01 


CH3OH 5 k \ V, 


= 


241.85 


BCD 


2.58 X 1.35 


80 


0.8 


0.03 


H2CS 7l,6 — > 6i_5 




243.98 


AB 


0.77 X 0.55 


8 


1.5 


0.01 


SO2 14o,l4 ^ 13i,,3 




244.25 


AB 


0.77 X 0.55 


8 


1.5 


0.01 


CH3OH 7t ^ 6t V, 


= 0, 1 


338.72 


comp.-ext. 


1.93 X 1.14 


86 


1.0 


0.1 



about 16.3° from the phase centre. The zenith opacities, mea- 
sured with the NRAO tipping radiometer located at the Caltech 
Submillimeter Observatory, were good during both tracks with 
t(348GHz) ~ 0.18 (scaled from the 225 GHz measurement 
via t(348GHz) ~ 2.8 x t(225GHz)). The receiver operated 
in a double-side band mode with an intermediate frequency 
of 4-6 GHz so that the upper and lower side band were sepa- 
rated by 10 GHz. The correlator had a bandwidth of 2 GHz and 
the channel separation was 0.8 125 MHz. Measured double-side 
band system temperatures corrected to the top of the atmosphere 
were between 150 and 500 K, mainly depending on the eleva- 
tion of the source. Details on the observational setup are given 
in Table |2] 

The initial flagging and calibration was done with the IDL 
superset MIR originally developed for the Owens Valley Radio 
Observatory and adapted for the SMAQ. T he imaging and data 
analysis was conducted in MIRIAD (Sault et al] Il995h and 
MAPPING. During the observations on November 11th, 2004, 
the position of the primary calibrator 3C111 was wrong in the 
catalogue by ~0.6". Therefore, we self-calibrated our secondary 
phase calibrator 0552+398 (16.3° from the source) shifting it in 
the map to the correct position. The solutions were then applied 
to IRAS 05358+3543. 



2.3. Single-dish observations witli tine IRAM 30 m teiescope 

In addition to the high resolution data, we mapped an area of 
70" X 180" in t he 5^^ ^ 4^ y, = CH3OH band, with the 
HERA receiver dSchusteret al.ll2004l) at the IRAM 30 m tele- 
scope in on-the-fly mode. The observations were performed in 
service-mode in February 2005, under excellent weather condi- 
tions (0.6-1.1 mm precipitable water vapour). The pointing was 
checked on Saturn and on a nearby source (0439+360) and was 
found to be accurate to ~ 6". Conversion from antenna tempera- 
ture to main-beam brightness temperature was performed by us- 
ing a beam efficiency of 0.4^ The beam of the 30 m telescope at 
241.8 GHz is ~ 10.2". Unfortunately, the overlapping in the UV 
plane of our PdBI and 30 m data is poor, since the PdBI base- 
lines start only only at 20 m, and we were not able to combine 
the two datasets and recover the short spacing information. The 
IRAM 30 m data are, however, used in the following discussion 
to study the extended structure of the methanol emission. 



' The MIR cookbook by Charlie Qi can be found at 
|http ://cfa-www. harvard. edu/~ cqi/mir cook.html | 
http://www.iram.fr/IRAMES/ 



3. Observational results 

In this section we present the results of the line observations per- 
formed towards IRAS 05358+3543. Since the observations were 
aimed at observing several bands of methanol transitions, the 
majority of detected lines comes from this molecule. However, 
emission from other molecular species with transitions close in 
frequency space to methanol is also detected. All detected lines 
are listed in Tables[3](for CH3OH) and|4](for the other molecular 
species); the dust condensations where the lines are detected are 
also reported in the tables. 

The large structured emission is sampled by the single dish 
methanol data, while the interferometric datasets allow us to 
zoom in on the gas around the protostars. These data range in 
angular resolution from 4" to 0.6", corresponding to 7200 and 
1100 AU, respectively. Another property of these data is that a 
broad range of structures are filtered out from the observations, 
due to the missing short-spacing flux, from ~ 25" in the BCD 
configuration of the PdBI at 3 mm to ~ 6" in its AB configura- 
tion at 1.3 mm. This means that the various datasets are differ- 
ently affected by the problem, and the comparison of transitions 
from different configurations is not straightforward. Moreover, 
the UV coverage is poorly sampled also at intermediate scales, 
and smaller structures can be filtered out. Comparing the single 
dish data for the 5^ ^ 4^^ v, = band of CH3OH to the inter- 
ferometric observation confirms that the PdBI data are missing 
fluxes. Although we cannot perform the same comparison for the 
rest of the dataset, negative features due to missing short spac- 
ings, heavily corrupt the other low energy line datacubes. Thus, 
the interpretation of the line interferometric data is potentially 
affected by missing flux problems. 

Although our observations are among the highest spatial res- 
olution studies of high-mass star forming regions existing to- 
day, there are obviously still limitations to what they resolve. We 
know from other observations that multiplicity on smalle r scales 
happens in star form ing regions of all masses (e.g., Weige ltet al.l 
1 19991: iMe^eath et ar.ll2005l) . The highest linear resolution of our 
observations (~ 1 100 AU) is reached only in one dataset, while 
the typical resolution of the other data is poorer. Within our own 
data, the two sources with the smallest separation, mmlb and 
mmla, are spatially resolved only in one dataset. Therefore, we 
cannot exclude that higher resolution would reveal more objects, 
and that the parameters we derive do not sample the gas around 
a single protostar. However, if this was the case, the multiple un- 
resolved sources would form from the same reservoir of dust and 
gas. Assuming that they are gravitationally bound, the parame- 
ters derived in our current study would still be relevant to the 
analysis of the collapsing cores. 
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R.A. [J2000] RA [J2000] 



Fig. 1. a:) In grey scale, the map of the integrated intensity of the 5o — > 4o v, = CH3OH-A line with the 30 m IRAM telescope 
(v = [-17, -16] km s levels from 0.5 Jy beam ' km s" in steps of 0.5). The dashed contours show the blue-shifted emission (v = 
[-23, -20] km s"'; levels from 0.5 Jy beam"' km s"' in steps of 0.5); the solid lines the red-shifted emission (v - [-14, -12] km s"'; 
levels from 1 Jy beam ' km s ' in steps of 0.5). The beam is indicated in the top left comer. The d otted circles outline the observed 
mosaic of seven fields. The pentagon locates the position of one of the H'^CO^ peaks observed bv lBeuther et al.l j2002al) . which is 
the candidate centre for the outflow traced by CH3OH. The soUd lines outline the region shown in panel b). b:) In grey scale (and 
grey contours), the map of the integrated intensity of the 5o — > 4o v, = CH3OH-A line with the Plateau de Bure Interferometer 
(v - [-24, -12] km s '). Level contours are from 1 Jy beam ' km s in step of 0.5. Dotted contours show the negative emission 
(-0.8 Jy beam ' km s '). The black contours outline the high velocity outflow (CO 2 — > 1 SMA data, Beuther priv. comm. Blue- 
shifted emission v - [-44, -24] km s"'; red-shifted emission v = [-8, 2] km s"'; levels from 30 to 100 Jy beam"' km s"' in steps of 
20). The four triangles mark th e positions of the main mm dust condensations. For illustration, the direction of the highly collimated 
outflow jBeuther et al ] l2002ai is shown by the arrows. For both flows, dashed lines are for the blue-shifted emission, solid lines for 
the red-shifted emission. The beam is indicated in the bottom left corner. 



3.1. Extended emission 

Emission from relatively low excitation lines is extended, as 
shown in Fig. 1(a), where the integrated intensity of the 5o — > 
4o, V, = CH3OH-A transition taken with the 30 m telescope is 
presented. 

For methanol, the emission at the cloud velocity shows a 
distribution very similar to the continuum emission. The off- 
set of the line emission peak from the continuum peak is real, 
since it was already observed in methanol and in other molec- 
ular species ( Beuther et al. 2002a). In their single-dish observa- 
tions, IBeuthgLetalJ (12002a|) detected red and blue non-Gaussian 
emission in the 5^: ^ 4^ methanol band, and suggested that both 
emissions are associated with a third outflow in the region to 
the west of the continuum peak. Their results are confirmed by 
our data, which show a bipolar distribution in the low excitation 
transitions of CH3OH (Fig. 1(a)) with a better signal-to-noise ra- 
tio. Within the pointing uncertainties and the low resolution of 



the data, the centre of this outflow seems to be associated with 
one of the H'-^CO^ peaks (indicated by a pentagon in Fig. 1(a)) 
detected bv lBeuflier et al.l (|2002a). If this were the case, our ob- 
servations would solve the problem of the powering so urce of 
the thi rd outflow of the region (C in the notation of Beuth er et al.l 
I2002al) . which was previously not assigned. However, higher an- 
gular resolution is needed to assign the centre of this flow beyond 
a doubt. 

The interpretation of the interferometric data of the same 
transition is not as straightforward as for the single-dish data. 
In Fig. 1(b), the integrated intensity of the 5o — > 4o CH3OH- 
A line is shown as seen in the PdBI observations, together with 
the direction of the two outflows originating from the vicinity of 
mml. The region west of the dust cores, where the bipolar dis- 
tribution is detected in the single-dish data, is only partially in 
the field of view of the PdBI observations. The intensity map is 
strongly affected by missing flux problems, which results in fil- 
tering out the extended emission. For this reason, the PdBI map 
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Table 3. Table of CH3OH detected transitions. In Col. 4, ★ is 
used to indicate the lines with extended emission; (?) for tenta- 
tive detections. 

Transition Rest frequency -Euppei Detections" 

(GHz) (K) 



5 

Table 4. Table of observed transitions from other molecular 
species, n Col. 4, ★ is used to indicate the lines with extended 
emission. 



Transition 


Rest frequency 


f 

^ upper 


Detections'' 




(GHz) 


(K) 




C b Z — » 1 


yo.4i j 


1 


1,Z,3,* 


L. b J — > 4 


o/i 1 n 1 A 
Z41.Ulu 






Crl3L.JN Ijg — > IZg 


zJK.y /Z 


ii 1 




L.rl3L.J>J 1^5 — > IZ5 


Zoy.Uzz 


1 ^o 




r^Tj 1 Q ^11 
L.rl3L.IN 1j4 — > IZ4 


z^y.UOJ 






L,rl3L.IN 1j3 — > IZ3 


ziy.uyo 


143 


l,line 


L.rl3L.iN LJ2 ^^2 


OQO 1 1 o 


1 no 




/^U /^"M 1 Q ^ 1 O 

L,rl3L.iN Iji — > 1Z[ 


z^y. 133 


/ 


l,line 


/''II r^xT 1 o ,10 


Tin 1 TO 
z3y. 13o 




l,line 


JrlLXJ Zi I — ^ Zi 2 


T/i 1 1 

Z41.J01 


0^ 






Z41.olu 


z4 


l,line 


xllNUU 112,10 ~^ ^^2,9 


1A 1 '70/1 
Z4 1 . / U4 






HNCO 112,9 — * 102,8 


241.708 


243 




HNCOllo,!! ^ 10o,io 


241.774 


70 




H2CS7i,6-^6i.5 


244.048 


60 




SO2 14o.i4 ^ 13l,l3 


244.254 


94 




C^'^S 7-^6 


337.396 


65 


l,line 


H2CS lOi 10 9i.9 


338.083 


102 




SO2 I84.14 ^ I83.15 


338.306 


197 




SO2 20u9 192.18 


338.612^ 


199 





" 1,2,3 indicate mml, mm2, and mm3, respectively. 1 is used for 
data where mm la and mm lb are not resolved; line is used for the 
position (1.2",0.6") off of mm2. 

* blend with the 7i ^ 61 v, = CHjOH-^ line at 338.615 GHz 



of this transition looks clumpy. The emission is detected all over 
the dust continuum peaks, often red- and/or blue-shifted. These 
emission spots could be caused by the high velocity CO outflow 
originating from the vicinity of mml (see Fig. 1(b)). However, 
we refrain from a further discussion of these features, given the 
missing short spacing information of our methanol data. 

The distribution of C^'^S line is puzzling (Fig.|2]l. While the 
2 — » 1 line extends over an area of approximately 20" x 20" 
around the main dust condensations, and the 7 — » 6 transition is 
found in association with mml (although marginally detected) 
and towards mm2, the 5^4 line is detected only on mml a 
and mmlb (see Fig. |2]i. In the 5 — > 4 C^''S line, two velocity 
components (see discussion on mml in ^3.2| l are detected, one 
at v/5r ~ -17.6 km s"' and the other at vi^r ~ -13 km s"', 
the first on mm la and the second on mm la and on the vicinity 
of mmlb. The emission of the 5^4 C^'^S line extends to the 
south-west, as does that of the 1 —> 6 line. However, the spectra 
are noisy, and the line is detected with a significance level not 
better than 2 cr. In Fig. |2] the level contours start from a value 



Cr Zl ' i ] Vl 
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of 3 cr, and the feature on the south-west of mm la is not seen in 
the 5 — > 4 transition. 

lOl mi & Cesaronil il999|) studied the excitation of CS and 
C^'^S towards a sample of ultra-compact Hii regions, and found, 
as expected given the energies of the levels involved, that the 
5 — > 4 C^^S transition is usually stronger than the 7 — > 6 line, 
making the non detection of the 5 — > 4 line on mm2 difficult to 
interpret. However, this line was observed in the AB configura- 
tions of the PdB interferometer, where source structures larger 
than 6" are filtered out. Therefore, the unexpected spatial distri- 
bution of this transition, as well as its profile, can be an artifact 
due to the missing short spacing information. 

Comparing the spatial distribution of the 2 — > 1 and 7 — > 6 
C"''*S lines, which are affected by missing flux problems on larger 
scales (>30" and >16", respectively), it emerges that they do not 
peak on the protostellar cores, but in the gas between the cores, 
where the temperature most li kely drops down. O bservations 
of the hot core in G29.96-0.02 tBeuther et alj|2007bi) . with the 
same frequency setup of our SMA data, show a similar ef- 
fect, with the 7 — > 6 C^^S line peaked in the gas between the 
(sub)mm continuum peaks. On the other hand, in G29.96-0.02, 
the Sg 1-1 ^'^SiO line peaks on the hot core. Since CS should 
desorb from dust grains relatively early in the evolution of the 
source dViti et al.l |2004 . while the SO and SO 2 abundances 
should increase with time, iBeuther et alj ('2007b) interpret the 
distributions of C^'^S and ^^SO in G29.96-0.02 as an evolution- 
ary indicator. In our case, we did not detect the Sg li ^"^SO 
transition toward IRAS 0535-1-3543. However, we observed four 
transitions of SO2, and all peak on the main dust condensation 
mml (Fig. [3]l. Although the evolutionary scenario outlined for 
G29.96-0.02 seems to be confirmed by our observations of the 
2 — > 1 and 7 — > 6 C-'^^S, and of the SO2 transitions, the distribu- 
tion of the 5 — > 4 C^'^S line does not find an easy interpretation 
within this theory. 

Similarly to C^"^S, also H2CS shows different morpholo- 
gies in the two transitions in our datasets (Fig. O, with the 
lOi.io — > 9i_9 line {Eup ~ 100 K) having a second peak at an 
offset position from mm2 (see the discussion below), while the 
7i,6 — > 61,5 is found only on mm la, despite its lower level energy 
{Eup ~ 60 K). We believe that these data are heavily contami- 
nated by missing flux problems at intermediate and short scales, 
and refrain from any analysis of this molecular species. 



3.2. Molecular emission from the dust cores 

The molecular spectra of the three main dust condensations, 
mml, mm2 and mm3, differ significantly from each other, re- 
flecting the different physical conditions, and probably time evo- 
lution, of the dust cores. Maps of the integrated intensity of 
several transitions are shown in Figs. [3] and |4] All high exci- 
tation lines peak on mml, and specifically on mml a for the 
dataset where the cores mmla and mmlb are spatially resolved. 
However, several moderately excited (Eup < 130 K) methanol 
lines are extended towards the second dust core mm2, and show 
a peak of intensity at (1.2", 0.6") from mm2. 

Source mml: Figure|5]presents the observed spectral bands 
at 238, 241 and 338 GHz, respectively, towards the main dust 
condensation mml. High excitation (E„p > 200 K) transitions 
are detected towards this core; however, the datasets where the 
two dust condensations mmla and mmlb are resolved, show 
that the emission from high excited lines is associated only with 
mmla. Only the 5 — > 4 C-^'^S transition (v ~ -13.4 km s"') is 
detected on mmlb. The lack of molecular emission from mmlb 
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Fig. 2. In grey scale the 1.3 mm continuum emission (from the 
AB configuration of the Plateau de Bure Interferometer). The 
black contours show the integrated line intensity of the 2 — » 1 
(a), the 5 ^ 4 (b), and of the 7 -> 6 C^'^S lines (c). Level 
contours are from 0.1 Jy beam"' km s"' in steps of 0.2 for 
the 2-^1 transition; from 0.3 Jy beam"' km s"' in steps of 
0.3 for the 5^4 fine, and from 1.5 in steps of 0.5 for the 
7 — > 6. Level contours for the continuum are in step of 0.003 
from 0.003 Jy beam"'. The dotted contours show the negative 
emission (-0.1 Jy beam"' km s"' for the 2^1 and 5 — > 4 lines; 
from -1.5 in step of 0.5 Jy beam"' km s"' for the 7 — » 6 tran- 
sition). The dashed lines in panel a) outline the region shown in 
the other two panels. The different beams are also given. 
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Fig. 3. Compilation of integrated intensity maps of transitions from several molecular species, shown as black contours. In grey 
scale (and in the solid contours in the first panel) the 1.3 mm continuum emission from the AB configuration of the Plateau de 
Bure Interferometer. Level contours for the line images are from 0.5 Jy beam"' km s"', in step of 0.3; for the 7i,6 — > 61,5 H2CS 
transition, from 0.3 Jy beam ' km s ', and from 1 Jy beam ' km s ' for the 184,14 — > 183,15 SO2 line. Level contours for the 
13o,i 12o,i CH3CN transitions are in step of 0.6 Jy beam ' km s '. The dotted contours show the negative emission: for the 
7i,6 — > 61,5 H2CS transition -0.1 Jy beam ' km s for the lOi io — > 9i,9 H2CS line -1.5 and -1.2 Jy beam ' km s"'. Continuum 
contours are the same as in Fig.|2l The beam is shown in the right bottom corner. 



could be a bias of the observations: the two dusts condensations 
are resolved only in the AB configuration of the PdBI, which was 
aimed at the detection of torsionally excited lines of methanol 
and had the 5 — > 4 C^'^S as the only low excitation transition. 
The two condensations have very similar properties for the con- 
tinuum emission (Beutheret al. 2007a}. However, other obser- 
vations suggest that mm la and mml b are indeed of different na- 
ture, as methano l class II masers (.Minier et alj2000l) . and a mid- 
infrared source dLongmore et alj 2006h are detected on mm la, 
but not on mmlb. A similar chemical differentiation of two cores 
of very close masses was recently found by Zhang et al. (2007) 
in the high-mass (proto)stellar cluster AFGL5142. As suggested 
for AFGL5142, a possible interpretation of our observations is 
that mmlb is in an earlier evolutionary phase than mmla, still 
not characterised by a rich chemistry. 

Source mmla: Several methanol transitions, in the ground 
state as well as in the first torsionally excited level, show two ve- 
locity components. The same profile is found in the 12^ lit 
CH3CN band (SMA data, Beutherpriv. comm.) and in the 13a — > 
12^ Hnes, in the 7i.6 61,5 H2CS line, and in the 5 ^ 4 C^'*S 
transition. HNCO is also showing the same behaviour in the the 

llo.ii — ^ 10o,i() line at ~ -17.8 and 15.3 km s"'. Two peaks 

are also detected in the 11 2,9 — > 102,8 transition, but the blend 
with the 5o 4o CH30H-E' line and with the ll2,io -> 102,9 



transition complicates the interpretation of this profile. On the 
other hand, the SO2 and HDO lines at ~ 241.6 GHz show only 
one velocity component, peaked at -16.6 and -14.9 km s re- 
spectively. For the HDO line, this corresponds to the velocity 
in between the two peaks detected in the other transitions. An 
example of the double peak profile of mmla is given in Fig. |6] 
where part of the 5ic — » 4^, Vt = 1 CH3OH band is shown. For 
the 7i,6 61,5 H2CS and the 5 ^ 4 C^'^S lines, this profile 
could be caused by missing flux; self absorption or missing flux 
problems could affect the low energy transitions of methanol. 
The hypothesis of self absorption is strengthened by the profile 
of HDO transition. However, since the double-peaked profile is 
found also at high energies, where the optical depth is unlikely to 
be high, even in rare molecular species, we believe that the two 
peaks are, at least for the torsionally excited lines of CH3OH, for 
CH3CN and HNCO, real, and not due to self absorption. In sev- 
eral cases, overlap between different lines complicates the iden- 
tification of the two components; however, a few lines (5o — > 4o 
and 5-1 4-i v, = 1 CH3OH; I34 Ma CH3CN) do not over- 
lap with any other transitions. By Gauss fitting their line profiles, 
we derive linewidths of ~3 and ~4. km s"', and system veloci- 
ties of ~ -17.6 km s ' and ~ -11.6 km s respectively. Both 
velocity components peak on mmla. 
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Fig. 4. Compilation of integrated intensity maps of CH3OH, 
shown as black contours. In grey scale (and in the solid contours 
in the first panel) the 1.3 mm continuum emission from the AB 
configuration of the Plateau de Bure Interferometer. Level con- 
tours for the line images are from: 0.3 Jy beam"' km s"' in step 
of 0.4, 1 Jy beam"' km s"' in step of 1, 1.5 Jy beam"' km s"' 
in step of 1, 0.2 Jy beam"' km s"' in step of 0.2, respectively. 
Dotted contours are used for negative values: for the Si — > 4-2-A 
line -0.2 Jy beam"' km s"'; -2 and -1 Jy beam"' km s"' for the 
7o — > 6q-E line; -2 Jy beam"' km s"' for the 7 ±2 — > 6±i-E line; 
-0.1 for the 5i ^ 4i CH30H-£' v, = 1 transition. Continuum 
contours are the same as in Fig. |2] The beam is shown in the 
right bottom corner. 
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Fig. 6. Spectrum of mmla at 241.2 GHz. The upper labels in- 
dicate the rest frequencies of the CH3OH lines for a velocity of 
-17.6 km s"', while the lower labels show the same transitions 
for a velocity of -1 1 .6 km s"' . 



We fitted the peak position of each velocity channels of the 
5-1 — > 4_i V, - I CHsOH-ii line, and detected a velocity gradi- 
ent along an elongated structure (Fig. |7(a)) . This linear structure 
seems to be shifted from the millimetre continuum peak. The 
same velocity gradient of the v, = 1 methanol lines is found in 
the 19 18 OCS line at 231 GHz (detected with the SMA, 
Beuther priv. comm.), and in the 6.7 GHz methanol maser tran- 
sition (Minier et al. 2000). In Fig. |7(a)[ we show for comparison 
the positions of the velocity channels of the maser line, and in 
Fig. 7(b)| the first moment map of the 19 — > 18 OCS transition. 
Moreover, observations with the SMA with an angular resolu- 
tion of ~ 3" (Beuther, priv. comm) show that the 2 1 C'^O 
line is a ssociated with the high velocity outflow ofl Beuther et al.l 
(I2002al) . and has a velocity gradient perpendicular to the linear 
structure detected in methanol. While an unresolved expanding 
or contracting shell of gas would also produce the double-peaked 
profile detected towards mmla, the linear velocity gradient de- 
tected in methanol and OCS cannot be explained by these sce- 
narios. 

The line profile detected towards mmla is reminiscent 
of the double-peaked profiles found in low and intermediate 
mass til ted circumstellar di sks, in optically thin and thick lines 
(iBeckwith & Sargenill993b . Moreover, the double-peaked pro- 
file arises from a region {R < 550 AU) comparable in size 
to the candidate disks surrounding high-mass (proto)stars (for 
a compilation of them, see Cesaroni et al. 2007). The detec- 
tion of a velocity gradient is not conclusive of a circumstel- 
lar disk. However, we believe that the velocity structure de- 
tected in methanol and perpendicular to the high-velocity out- 
flow, together with the double-peaked profile, and the size of 
the emitting gas are strongly suggestive of a rotating struc- 
ture in IRAS 05358-mmla. To investigate whether this struc- 
tur e is a circumstellar d isk or a toroid like the ones discussed 
bv lCesaroni et al.l (|2007|) . higher linear and velocity resolutions 
are needed to study the velocity field of the innermost gas around 
mmla. We note however that toroids have typical radii of several 
thousands AU, while the radius of the structure around mmla 
has an upper limit of 550 AU. 

Alternatively, two different hot components with a separa- 
tion of less than 1 100 AU would also produce a double-peaked 
profile, and a velocity gradient along a linear structure. Since the 
two lines have very similar intensities and linewidths, the physi- 
cal conditions in the two cores would be comparable (as derived 
in paragra ph §4.jj i^The total mass of the binaries can be derived 



ra p 

via (e.g!. IChen"iraLll2006l) 



M,„, = 



A X Av^ 
G X sirfi{i) 



(1) 



where A is the semi-major axis of the elUptical orbit of the two 
stars, G the gravitational constant, / the inclination angle, which 
is unknown, and Av the relative velocity between the two cores. 
An upper limit to the major axis is given by the resolution of the 
data (0.6", corresponding to ~ 1 100 AU, at this distance). In sec- 
tion ^4.11 the sizes of the two sources are derived by fitting the 
13t ^ 12^ CH3CN band. Therefore, the lower hmit to the major 
axis of the orbit is taken as the sum of their radii, correspond- 
ing to 0.2" (360 AU), which is in agreement with the extension 
of the linear structure in Fig. |7(a)| Hence, the total mass of the 
system ranges between 22 and 7 Mo/sin^/, in good agreement 
with the value of 1 3 M p derived from the Lyman continuum flux 
dBeuflier et alj|2007 a'). 

In both cases, a rotating structure or two close-by hot cores, 
the single-peaked line profile of HDO would suggest that the 
emission is coming from a larger envelope of gas. 
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Fig. 7. a: Presented is the position-position (P-P) diagram of the 
5-1 — > 4_i V, - 1-E CH3OH line. The star marks the position 
of the 1.3 mm continuum peak. The values above each mark 
are the LSR velocities of the corresponding channels. The error 
bars represent statistical errors of the fits. The error on the posi- 
tion of the continuum peak are negligible. The black dots on the 
east mark the po sitions of the veloc ity channels of the 6.7 GHz 
methanol maser dMinier et al The ofl'set between the two 

linear structures can be accounted for pointing uncertainties, and 
error in the absolute position of the VLBI data, b: First moment 
map of the 19 ^ 18 OCS transition. The white triangles mark 
the positions of the dust condensations mmla and mmlb. The 
beam is shown in the bottom right corner. 



To discriminate between the different possibilities discussed 
above, higher spatial resolution observations will be required to 
resolve and image the system in more details. Among the cur- 
rent cm-mm interferometers, only the VLA in its A configuration 
would allow observations with a spatial resolution significantly 
higher than the one of our current data. Otherwise, studies of 
such systems at resolutions of 0.1" or less in spectral lines other 
than NH3 will be possible in the submm windows only when the 
Atacama Large Millimeter Array (ALMA) will come on line in 
the next few years. 

Source mm2: the molecular spectrum of mm2 shows 
emission from several species, with moderately excited 
{Eup < 130 K) methanol lines. Linewidths are compara- 
ble to the one of the components in mml, and emission 
is at v=-15.3 km s^'. However, stronger emission is de- 
tected in several methanol lines at (1.2", 0.6") off of mm2 
(q'2ooo=05''39"'12'.86, 52()(K)=+35°45'51".9), at a position where 
no continuum emission is detected (see Fig.|4]l. At the same po- 
sition, the CH3CN k = 0,1,3 lines are also detected, while the 
7 — > 6 C"^^S and lOi io 9i,9 H2CS transitions have a secondary 
emission peak in the vicinity of this position. Linewidths are of 
the order of 3.7 km s"', at the same velocity of mm2. 

To rule out the hypothesis that the diff'erent distribution of 
the gas with respect to the dust is not real, but due to optical 
depth effects, the intensity of the lines should be mapped in the 
wings, where the optical depth is lower. Unfortunately, the sig- 
nal to noise in the CH3CN spectrum and in the SMA data, which 
are the two datasets less affected by line-overlapping problems, 
is not good enough to perform this check. However, given the 
range of energies of the lines showing a second peak between 
mml and mm2 (see Fig [3] and HJ, we believe that this spatial dis- 
tribution is real. In Fig. [8] the CH3CN and CH3OH spectra are 
shown towards the offset position. In the rest of the paper, we 
will refer to this position as to mm2-line. The molecular spec- 
trum of mm2 looks very similar to the one of mm2-line, but line 
intensities are usually a factor of 1.2 weaker than on mm2-line. 
Since the two positions have a separation of roughly half the 
beam size of our observations, we cannot rule out that part of 
the emission at mm2 is coming from mm2-line. 

A possible interpretation for the molecular emission at mm2- 
line is that it is caused by sputtering or thermal evaporation of 
the grain mantles, followed by gas-phase reactions, due to the 
interaction with one or more of the outflows in the region (see 
Fig. |9j. Although mm2-line is characterised by narrow features 
{Avisfi ~ 3.2 km s"'), the emission is red-shifted respect to 
mml, thus strengthening the outflow scenario. However, the de- 
tection of a high-density tracing molecule like CN3CN at this 
position does not easily fit with this hypothesis. Alternatively, a 
low-mass dense core could be responsible for this emission. In 
this case, an upper limit to the mass of mm2-hne can be com- 
puted from the continuum emission at 1.3 mm, assuming a tem- 
perature of 60 K as derived from the analysis of its methanol 
spectrum (see § 14.31 1. This results in a mass of Me, ~ 0.08 Mq. 
In § 14.31 we derive a lower limit of 10"^ for the abundance of 
methanol at mm2-line, whi ch is comparable to the values found 
in low-mass star hot cores (iMaret et ani2005h . These cores are 
however usually warmer than 60 K. Otherise, the enhancement 
of the CH3OH abundance can also be explained within the hy- 
pothesis that the emission is caused by the interaction of molec- 
ular outflows with the ambient molecular cloud. 

Source mm3: the core mm3 presents the simplest molecu- 
lar spectrum among the several dust condensations in the region. 
The emission associated with mm3 is extended and not com- 
ing from a compact core, as all the lines detected are seen as 
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Fig. 9. In grey scale, the integrated intensity of the 
72 ^ 62 Vt = CH30H-£ with the SMA (from 
2.5 Jy beam ' km s ' in step of 1), overlaid on the CO 
2 — > 1 high velocity outflow (SMA, Beuther, priv. comm. 
Blue: solid lines v - [-44,-24] km s"'; red: dashed lines 
V = [-8, 2] km s"'; levels from 1.5 Jy beam"' km s"' in steps of 
1). The four triangles mark the positions of the main mm dust 
condensations; the square outUnes the mm2-line position. 



negative features, an artifact due to missing spacing s. Therefore, 
only a q ualitative analysis of these data is possible. Beuther et alj 
( l2007al) computed its spectral energy distribution and derived 
an upper limit to the dust temperature of ~ 20 K, a value that 
is constrained by the weakness of the continuum emission at 
690 GHz. In paragraph 94.31 we discuss the molecular spec- 
trum at this position, and conclude that the gas around mm3 is 
cold (< 30 - 40 K), and dense (10^ - - 10'' cm"^). Therefore , our 
results agree with the conclusion of i Beuther et alj (12007 a') that 
mm3 is a candidate of a cold massive core in an early evolution- 
ary stage. Our observations of CH3OH and C^^^S suggest that 
these molecules have an extended, flat distribution, and show 
no trace of the structure seen in the dust continuum emission, 
a result reminiscent of observations of low-mass starless cores 
(e.g.,|Caselli et al. 1999). Even H' ^CO^, a good tracer o f quies- 
cent gas, shows no peak on mm3 dBeuther et alj|200^ah . These 
results suggest that mm3 is indeed a candidate massive starless 
core. However, high angular resolution observations that trace 
the dense gas and investigate potential signs of star formation 
activity are needed to test this hypothesis. 



4. Derivation of physical parameters 

4.1. CHiCN 

The physical properties of the gas around mml are derived by 
the analysis of the 13^ —> \2k CH3CN lines. The analysis was 
carrie d out with the Xclass program (discussed in lComito et alj 
l2005h . which uses an LTE model to produce synthetic spec- 
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Fig. 10. Spectrum of the 13^: 12^ CH3CN band towards the 
main dust condensation mml; overlaid on the data, in black, 
the best fit synthetic spectrum. All CH3CN lines are labelled in 
the upper axis for the velocity component v - -17.6 km s 
The frequencies of the corresponding CH^^CN transitions are 
labelled in the lower axis. 



tra, and compares them to the observations. The molecular data 
are from the CDMS (Mufler et al. 2001) and JPL (Pickett et al] 
1998) databases. The parameters defining the synthetic spectrum 
are: source size, rotation temperature, column density, velocity 
width and velocity offset (with respect to the systemic velocity 
of the object). Several velocity components, which are supposed 
to be non-interacting (i.e. the intensities add up linearly), can be 
used. 

Given the low signal-to-noise ratio of the spectra, it is possi- 
ble to perform this analysis only on mml. For this position, we 
assumed that all the CH3CN transitions trace the same warm, 
dense gas around the central object and come from the two ve- 
locity components detected in several molecular species (see 
paragraph 93. 2t , although the lower excitation lines may have 
contribution from a more extended component (Fig. 13). Since 
the line profiles consist of two peaks of the same intensity, we 
used the same physical conditions for the two components. The 
source size is degenerate with temperature in the case of com- 
pletely optically thick lines, and with column density for com- 
pletely optically thin lines. In our case, the ^ > 3 lines are 
optically thin, while the k < 1> are optically thick. Therefore, 
the degeneracy between source size and column density can be 
solved. Constrains on the column density of CH3CN come also 
from the non detection of the \3k 12jt CH^^CN band, which 
is very close in frequency to the CH3CN lines. We found that 
the source size is constrained to ~ 0.2" for each components, 
in agreement with the size of 0.2-0.3" of the linear structure of 
Fig. |7(a)| For this value, the best fit corresponds to a tempera- 
ture of 220 K, and a column density of ~ 4 x 10'^ cm"^ per 
velocity component. The uncertainties on these values are re- 
ported in Table |5] Assuming that the CH3CN emission comes 
from mml a, and using the hydrogen column density derived by 
[Beuther et al. (2007a) for the AB configuration data at 1.2 mm 
(A^Ht ~ 2.2 X 10^^ cm"^), this results in a total abundance of 
CH3 CN relative to H2 of 4 x 10 w hich is typical of hot c ores 
(e.g.. lHatchell et al.lll998l) . However. iBeuther et al.1 (l2007al) de- 
rived the H2 column density for a temperature of 50 K, while our 
estimate of Tk is higher. The abundance increases to 2 x 10"^, 
when re-computing the H2 column density of mm la at 220 K 
(A^H, ~ 4.6 X 10^^ cm-2). 
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4.2. SO2 

Four SO2 transitions were detected towards the main condensa- 
tion mml. They range in energy of the upper level from ~ 24 
to 200 K; therefore, although they look point-like in all observa- 
tions, they likely trace different gas around mml. We analysed 
the SO2 emission by means of the LTE technique discussed in 
paragraph 34.11 Since the lines do not show any trace of the 
double-peaked profile detected in other molecular species, we 
used a model with one velocity component for SO2. Assuming 
a source size of 1", the best fit is found for T - 150 K and 
N{S02) = 10'^ cm"^, which corresponds to an abundance of 
SO2 relative to molecular hydrogen of 2 x 10"*^ using Tdust - 
50 K (A^H, ~ 6.6 X 10^3 cm-2), and of ~ 5 x lO"*^ for a tem- 
perature of 150 K (Ah, ~ 2 X 10~^ cm"-). The errors on the 
parameters are listed in Table |5] 



4.3. CHiOH 

For the analysis of the CH3OH spectrum, we used a technique 
similar to the one discussed f or CH^CN 04.1ll and imple mented 
to the LVG approximation ("Leurin i et al.l |2004 |2007|) . Since 
methanol is detected at several positions, in the following discus- 
sion the properties of each dust condensation will be discussed 
separately. The errors on the parameters are listed in Table|5]for 
each condensation. 

Source mmla: In our analysis of mmla, we included only 
emission coming from the - 1 lines, as their optical depth is 
lower than for the ground state, and their emission is confined 
to the gas around the dust condensation, while the v, = transi- 
tions are more extended, and are affected by problems of missing 
flux. Moreover, the v, = 1 lines have similar level energies, while 
the ground state transitions in our data range from a few K to 
~ 200 K, thus not tracing the same gas. The model includes the 
three torsionally excited bands observed, the 2^ 1;; at 3 mm, 
the 5k — » 4i at 3 mm, and the Ik at 0.8 mm. Infrared 

pumping by thermal heatin g of the dust is in cluded in the form 
of an external grey-body (L eurini et al.ll2007l) . with yS = 1.6 and 
'rioo/jm=l-5. This is necessary to explain the excitation of the tor- 
sionally excited lines, as, with very high critical densities (10''^- 
10" cm"^) and high level energies (T > 200 K), they can be 
hardly populated by collisions, but trace the IR field instead. 
On the other hand, the ground state li nes can be populated by 
collisions and by the infrared pumping. iLeurini et al.l (|2007|) in- 
vestigated the effects of infrared pumping on the excitation of 
meth anol, and foun d that, as in the case of other molecules (e.g. 
CS, ICarroll & Go ldsmith 1981), the infrared pumping mimics 
the excitation by collisions of the v, = lines. This results in a 
degeneracy between the density and the infrared field of the ther- 
mal dust which makes any determination of the density through 
the analysis of the ground state lines of methanol impossible, for 
sources like hot cores. As in the discussion in paragraph 34.11 
we used two velocity components to model the data. 

The best fit is reached for values very similar to the ones 
derived with CH3CN, source sizes of 0.15" and a temperature 
of 220 K. As the line intensities are very similar in the two ve- 
locity components in almost all lines, the same parameters are 
derived for both of them. The column density for methanol is 
equal to ~ 2 x 10'^ cm"^ per component, which, using the H2 
column density derived by Beuther et al. (2007a) for mmla for 
50 K, corresponds to a total abundance of methanol relative to 
molecular hydrogen of 2 x lO"*", or to 10"^ if we adopt 220 K 



as dust temperature. Both v alues are typical of massive hot cores 
(iMenten et al.lfT986l fT988h . 

The fit well reproduces the intensities of the 7 = 5,7 series, 
but it fails to correctly predict the strength of the lines at 96 GHz, 
which are weaker in the model than in the observations. A plau- 
sible explanation is that 3 mm lines come from a more extended 
region than the others, and therefore the real beam dilution is 
less than that obtained in the model with a source size of 0.15". 
With the resolution of our data, we cannot verify this hypothesis. 
The 2k — > U lines were observed in both the BCD and the AB 
configurations of the array, and no flux is lost in the data with the 
higher resolution. However, the CLEAN beam of the AB data at 
3 mm is significantly larger than the one at 1 .3 mm, and different 
emission sizes between the 3 mm and the 1 .3 mm lines cannot be 
ruled out. Moreover, the 96 GHz lines peak at 300 K in energy 
of the upper level, while the other bands are at higher energies 
(330-530 K), and they are therefore more likely to trace the gas 
closer to the central heating source. 

Source mm2: Since no torsionally excited lines are detected 
toward mm2, the analysis on the methanol spectrum was per- 
formed on the ground state lines, and no external radiation field 
was used in the calculations. As discussed for mmla, the anal- 
ysis of these lines is not trivial, as they range in energy of the 
upper level from 35 K to more than 100 K (see Table |3j, and 
therefore trace different regions of gas. Moreover, some v, = 
lines have extended emission and are, thus, affected by miss- 
ing fluxes problems. The analysis of the datacubes at 241.7 and 
338 GHz shows that several methanol lines have a spatial distri- 
bution similar to that of the 5o 4() v, = CH3OH-A transition 
(Fig. 1(b)). We excluded these lines from our study of mm2, and 
limited the analysis to the transitions that are peaked on the main 
dust cores. An asterisk in Table |3]marks the lines with extended 
emission, excluded from the analysis of mm2. 

The continuum data (Beuther et al. 2007a) reveal a complex 
area around mm2, with several sub-sources which may not be 
of protostellar nature, but which could be caused by the out- 
flows in the region. Given the relatively low resolution of the 
observations of the methanol ground state lines, the molecular 
emission detected toward mm2 could arise not only from the gas 
associated with the dust core mm2a, but also from the other sub- 
sources around it. Therefore, in modelling the molecular spec- 
trum of mm2, we used the beam size as source size, and derived 
the average physical parameters of the gas over the beam. This 
approach is justified by the low optical depth of the methanol 
lines for the physical conditions found in the analysis. The same 
technique was used to model the methanol emission at the posi- 
tion (q'200o=05''39'"12'.86, 52(K)0=+35°45'51".9), north-west of 
mm2. Since no continuum emission is detected there, the nature 
of the molecular emission is unknown. The spectra toward the 
two positions, mm2 and mm2-Une, are very similar: emission 
comes from the same transitions, which have similar linewidths, 
and emit at similar velocities. However, the emission from mm2- 
line is stronger than on mm2 in all methanol lines, hence sug- 
gesting that the methanol abundance is enhanced at this position. 
The best fit is found for a temperature of 60 K, and a column 
density of 4 x 10^ cm"- and 7x10^ cm"^ for mm2 and mm2- 
line, respectively. The column densities at the two positions are 
4x10'^ cm"- for mm2, and 7x10'^ cm"^ for mm2-line. The cor- 
responding abundances of methanol relative to H2 can be calcu- 
lated by using the H2 column density derived by iBeuther et al.l 
(l2007al) from the BCD configuration data. The methanol abun- 
dance of mm2 corresponds to 6 x lO"'^, a value closer to the ones 
derived towards t he inner regions of low-mass protostars (e.g., 
iMaret et aLll2005h . than to values found in high-mass protostars 
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dMenten et alj|1986i [T988I) . In §[321 we noticed that, given the 
separation between mm2 and mm2-line, at least part of the emis- 
sion of mm2 could come from the offset position. Therefore, the 
CH3OH column density derived for mm2 should be regarded as 
an upper limit to its true value. 

An upper limit to H2 column density of mm2-line can be 
computed from the continuum emission at 1.3 mm, N ~ 5 x 
10^^ cm"^ at 60 K. This results in an abundance of 1 x 10"^ 
for methanol, which is a factor of 17 higher than on mm2. 
However, a different dust temperature would result in different 
values for Nh^ and Xch,oh- The enhancement of methanol is 
usually ascribed to evaporation of icy grain mantles in the vicin- 
ity of warm embedded sources, e.g. hot cores, or to desorption of 
grain mantles in sho cks associated to molecular outflows (e.g., 
lBachilleretal] |1995V Since no embedded sources are detected 
in the continuum emission at this position down to an upper limit 
of 0.08 Mq, it seems likely that the molecular emission associ- 
ated with mm2-line is related to the outflow activity of the re- 
gion. 

In the case of mm2, and mm2-line, the temperatures derived 
from m ethanol are comparable to those adopted by Be uther et al.l 
(l2007ah for the analysis of the continuum emission; therefore, 
only one value for the methanol abundance is given at this posi- 
tion. 

Source mm3: As listed in Table |3] only molecular transitions 
from the ground state, and with upper energies of less than 70 K 
are detected toward mm3. Since all the lines detected toward this 
core are seen in absorption due to missing fluxes, with emission 
only in the 5o — > 4-o-A and 5_i — > 4_i-£' transitions on top of the 
negative features, only a qualitative analysis of the data is possi- 
ble. We ran LVG models for several temperatures and densities, 
to find the range of parameters compatible with our observations. 
The detection of the 5+2 — > 4±2-£ lines, and the non detection 
of the higher k transitions constrain the H2 density to the range 
~ 10^-10^ cm"-', while the tentative detection of the k = Q-A 
and k = -1-E lines in the 7^ —> 64 band, but not of the other 
transitions, sets an upper limit to the temperature of 30-40 K. 



4.4. Uncertainties on tlie piiysical parameters of tlie gas 

Error ranges for the various parameters can be estimated by a^^ 
analysis, although this can only give errors within the assump- 
tions of the model. It cannot assess errors due to the LTE and 
LVG assumptions and to the treatment of the gas as consisting 
of a finite number of non-interacting components of homoge- 
neous conditions. As for the LTE assumption, densities in hot 
cores usually are in a range which makes it reasonable. 

The analysis of CH3CN and CH3OH infers a source size of 
~ 0.2" for the gas responsible of the emission of the two velocity 
components found toward mmla. This value is relatively well 
constrained by the CH3CN band, which shows optically thick 
and thin lines, and by the non detection of CH^^^CN. Since this 
result is in agreement with the size of the linear structure found 
from the fit of the positions of the velocity channels of methanol 
(see Fig. |7(a)| i, we fixed the source size to this value, and de- 
rived confidence levels for the estimates of the temperatures and 
column densities. However, one should remember that for op- 
tical thin lines, like the v, = 1 CH3OH transitions, the source 
size and column density are degenerate parameters, and the es- 
timates of the column densities would change, if the source size 
were different than our assumption. For CH3OH toward mmla, 
since these data were taken with a beam of ~ 0.7" and we as- 
sumed a size of 0.15" in the calculations, the CH3OH column 



density can only be a factor of 10 smaller than indicated by our 
results. For a source size of 0. 15", good fits to the data are found 
also for lower temperatures, but higher column densities. For 
T ~ 100 K the model requires column densities of > 10'^ cm"^, 
which would correspond to very high abundances relative to H2 
for such temperatures (Xch,oh ^ 10"^). For temperatures higher 
than 150 K, the methanol column density is confined to the in- 
terval 4xl0'^-2xl0'^ cm"^. A lower Hmit to the kinetic 
temperature is found at ~ 75 K, where the torsionally excited 
lines of methanol are not anymore efficiently pumped. 

Similarly, the column density of SO2 on mml, and of 
CH3OH on mm2 and mm2-line would increase if they were con- 
strained to a region smaller in size than the value used in the 
calculation. For example, the column density of SO2 would in- 
crease of one order of magnitude by using a sorce size of 0.2" 
as derived from the analysis of CH3CN and CH3OH. 

The 3 cr confidence levels derived for each parameter from 
the analysis are given in Table |5] together with their best fit 
values. 



5. Various stages of star formation 

The analysis of the spectral line emission in IRAS 05358 reveals 
four dust condensations with a maximum projected distance of 
18 000 AU. At least three of them show signs of active star for- 
mation. The core mmla appears to be the most evolved source of 
the region, associate d with a hyperc ompact Hii region, and class 
II methanol masers. iBeuther et al.l (|2007a) derived the Lyman 
continuum flux of mmla corresponding to a Bl ZAMS star of 
13 M0, and a large reservoir of accreting material (~ 10 M©). 
However, this estimate is obtained from relatively low resolu- 
tion data, where the cores mmla and mmlb are not spatially 
resolved. From high resolution data, which are however affected 
by missing fluxes, the estimate to the mass of mmla is 1 Mq, and 
0.6 Mq for mmlb. Our observations suggest that mmla also har- 
bours a candidate circumstellar disk. A rough estimate of the lu- 
minosity of mmla can be derived through the Stefan-Boltzmann 
formula, assuming a radius of 180 AU, and a temperature of 
220 K. By adding the values for the two components, the total 
luminosity ofmmlaisL ~ 6000 Lq, in agreement wit h the lumi- 
nosity of 10^^2 Lo of a Bl ZAMS star (lPanagiall973h . Since the 
total luminosity of the cluster is 6300 L©, this result supports the 
int erpretation t hat mm la is the main powerhouse of the region. 

IViti et all (I2OO4 derived the abundances of several 
molecules in the gas phase, as function of the age, and of the 
mass of the heating central object. The abundances of CH3OH, 
CH3CN, and SO2 of mmla are in agreement with a 15 star 
of approximately 10"^^ yr, although the best fit value for the col- 
umn density of CH3OH is higher than the theoretical values, but 
still consistent with the models within the uncertainties of the fit. 
Decreasing the mass of the heating source in the models results 
in increasing its age. The value of 10"*'^ yr is in very goo d agree- 
ment with the value obtained by iBeuther et al.l ( l2002al) for the 
dynamical timescale of the molecular outflows originating from 
the vicinity of mml (f ~ 37 000 yr). 

The physical parameters of the gas derived from our analy- 
sis, (temperature and abundances) are similar to hot cores, thus 
implying that mmla has already reached this evolutionary stage. 
However, its mole cular spectrum diffbrs from the one of typical 
massive hot cores. IBeuther e t aP ('2007b') studied the hot core in 
G29.96-0.02 (L ~ 9 x lO'* Lq, Olmi et al. 2003) with the SMA, 
with the same frequency setup we used for our SMA observa- 
tions. Although the linear resolution reached for G29.96 is com- 
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Table 5. Overview of the physical parameters in IRAS 05358+3543: - indicates parameters which cannot be derived with our data. 
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' derived from the 1.3 mm continuum data presented bv lBeuther et al 
* the values refer to the total column density of the two velocity components 



1 2007 al) and using the temperature derived from each molecular species 



parable to ours, the spectrum of G29.96 shows stronger methanol 
lines than IRAS 05358+3543 and higher excitation transitions; 
moreover, molecular species which are detected in G29.96 are 
not in our dataset. The same is found when comparing our obser- 
vations with Orion-KL (L ~ 10^ Lq), even with the po or angular 
resolution of the CSO telescope jSchilke et al.lll997h . Possible 
explanations for these differences are that mml a is indeed a) 
in an earlier evolutionary stage than the other two sources, at 
the beginning of the hot core phase, and will eventually reach 
their chemical richness; or b) hosts a less luminous hot core, 
which will never produce the rich chemistry of G29.96-0.02 
or Orion-KL. Both interpretations render mm la an interesting 
source, since they collocate it in a different region of the age- 
luminosity plane of hot cores, which still has not received much 
attention. Studies of large samples of hot cores at the same hnear 
resolution, and with the same frequency setups (similar to that of 
lBeutheij[2007l but on a larger sample of sources) are necessary 
to better understand this evolutionary phase of star formation. 
Particularly interesting is the comparison with observations of 
hot cores in the vicinity of low- and intermediate-mass proto- 
stars, which could clarify whether the poorness of the chemistry 
in IRAS 05358+3543 is due to its age or it is rather an intrinsic 
property of less luminous hot cores. 

Although the mm-wavelength thermal continuum emission 
of mmlb is similar in intensity and extent to that of mml a, even 
if not as strong, mmlb is not associated with any cm contin- 
uum emission down to a threshold of 1 mJy, while mml a is. 
Moreover, the molecular spectrum of mmlb seems very differ- 
ent from that of mml a. Of all the species we observed, only 
C^'*S is clearly found in mmlb. This is not a bias of our observa- 
tion as other observations confirm our results. Class II methanol 
masers and a mid-IR source are detected only towards mml a, 
hence confirming the different nature of the two cores. The core 
mmlb could be at an earlier stage of evolution, when the heating 
from the central star still has not produced any complex chem- 
istry, nor formed a hypercompact Hii region. 

The region around mm2 represents a challenge for the in - 
terpretation of the continuum emission (iBeuther et al.ll20()7ab . 
as well as for the spectral data. Several sub-sources are iden- 



tified around the main dust condensation, which are probably 
caused by the interaction of molecular outflows with the molec- 
ular cloud. Their total mass corresponds to 5 M©. The molecular 
emission associated with the dust core does not show any sign 
of highly excited lines, but points out at a warm (T ~ 60 K), 
dense (n > 6 x 10'' cm"^) gas with a methanol abundance of 
~ 7 X 10"^. The XcH,OH of mm2 is more similar to the values 
derived in star forming regions of lower mass than to high-mass 
young stellar objects. Given the total mass of the region, and the 
low abundance of methanol, we speculate that mm2 is a low- 
intermediate mass protostar. However, the spatial resolution of 
our data does not allow us to verify whether or not at least part 
of the emission from mm2 is indeed not coming from mm2-line. 

The coldest and least active core of the region is the dust con- 
densation mm3. Its mm continuum emission shows that mm3 
is a compact, unresolved so urce even at the hig hest resolution 
reached by our observations (iBeuther et al.ll2007 a). No compact 
molecular emission is detected from this core, but all the lines 
observed toward it are seen as negative features due to the fil- 
tering out of large structur es of the data. The analysis of the 
spectral energy distribution ( Beuther et al.l2007ah . and of the gas 
around the core confirm that mm3 is a dense (n> 10^ cm"^), 
cold (T < 40 K) condensation, surr ounded by a large res ervoir 
of accreting material (M ~ 19 Mn. IBeuther et aLll^OTah . The 
extended, flat distribution of molecules like CH3OH, C-'^S and 
H'^CO^, with no trace of the structure seen in the dust contin- 
uum emission, is reminiscent of low-mass starless cores, where 
C-bearing molecules are freeze-out onto the grain man tles in a 
cold and dense environments (e.g.. ICaselli et al.lll999h . Could 
mm3 be a massive starless core? Massive starless cores can be 
found in the vicinity of active sites of massive star-formation, 
and should be quiescent displaying no signposts of star forma- 
tion. Its properties suggest that mm3 is indeed a massive cold 
core in a very early evolutionary stage. However, with the obser- 
vations we currently have we cannot exclude that star formation 
have already started in this core. In particular, we cannot exclude 
that mm3 is the driving source of one of the outflows detected in 
the region. 
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6. Summary 

Our new interferometric data resolve at least four cores in 
the high-mass protocluster IRAS 05358-H3543. By analysing the 
molecular spectrum of each condensations, we characterised the 
properties of the gas surrounding it. Our main results are sum- 
marised in the following: 

- the main powerhouse of the region, mm la, harbours a hot 
core with T ~ 220 K, and the central heating source has a 
chemical timescale of lO'*-^ yr. Our data suggest that nrnila 
might host a massive circumstellar disk; 

- although the properties of the mm continuum emission of 
mmlb are very similar to the one of mmla, the two sources 
differ significantly in the cm and mid-infrared spectrum. This 
core could be in an earUer stage of star formation than mmla, 
since no molecular emission is detected toward it, with the 
only exception of the 5^4 C^^^S transition; 

- given the low abundance of methanol, mm2 could be a low- 
intermediate mass protostar; 

- strong emission is detected in several molecular species to 
the north-west of mm2, at a position where no continuum 
emission is detected. We suggest that this is caused by the 
interaction of the outflows with the ambient molecular cloud; 

- the least active source, mm3, could be a starless massive 
core, since it is cold (T < 20 K), with a large reservoir of 
accreting material (M ~ 19 Mg), but no molecular emission 
peaks on it. 
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